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Abstract 

The recent F2 data can be described using the DGLAP evolution equations 
with an appropriate choice of input distributions and the choice of the start- 
ing scale for the Q 2 evolution. We demonstrate in this paper that we cannot 
conclude that there are no significant shadowing corrections at HERA kine- 
matic region using F% data only. In this paper we calculate the shadowing 
corrections to the longitudinal structure function Fl(x, Q 2 ) and to the charm 
component of the proton structure function F%(x, Q 2 ) at the HERA kinematic 
region using an eikonal approach. We demonstrate that the shadowing cor- 
rections to these observables are very large and that the charm production is 

strongly modified at small- 2;. Our results agree with the recent few HI data. 
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I. INTRODUCTION 



The parton high density regime in deep inelastic scattering is one of the frontiers in 
perturbative QCD (pQCD). This corresponds to the small Bjorken x region and represents 
the challenge of studying the interface between perturbative and non-perturbative QCD, 
with the peculiar feature that this transition is taken in a kinematical region where the strong 
coupling constant a s is small. The transition between the perturbative regime and the non- 
perturbative regime may provide an understanding of the soft interactions in the language 
of QCD and establish whether the Regge model can be justified from QCD principles. From 
the experimental side, data from electron-proton collider HERA at DESY (-^/s = 300 GeV), 
provide new insights into the structure of the proton. The high energy allows to probe the 
proton in hitherto unexplored kinematical regions. The collider HERA allows the measure 
of structure functions F 2 , Fg, F L and the cross sections of vector meson production and 
diffractive processes for Bjorken x values down to 10~ 4 . 

In the region of moderate x (x > 10 -2 ) the well-established methods of operator ex- 
pansion and renormalization group equations have been applied successfully. The DGLAP 
equations [|TJ, which are based upon the sum of QCD ladder diagrams, are the evolution 
equations in this kinematical region. However, in the small x region, the density of gluons 
and quarks becomes very high and a new dynamical effect is expected to stop the further 
growth of the structure functions. In this kinematical region we are dealing with a system of 
partons which are still at small distances, where a s is still small, but the density of partons 
becomes so large that the usual methods of pQCD cannot be applied ||. About seventeen 
years ago, Gribov et al. || have performed a detailed study of this region. They argued 
that the physical processes of interaction and recombination of partons become important 
in the parton cascade at a large value of the parton density, and that these shadowing cor- 
rections could be expressed in a new evolution equation - the GLR equation. This equation 
considers the leading non-ladder contributions: the multi-ladder diagrams, denoted as fan 
diagrams. The main characteristics of this equation are: it predicts a saturation of the gluon 
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distribution at very small x; it predicts a critical line, separating the perturbative regime 
from the saturation regime; it is only valid in the border of this critical line. Therefore, 
the GLR equation predicts the limit of its validity. In the last decade, the solution |3]-|6| 
and possible generalizations [0,|J of the GLR equation have been studied in great detail. 
Recently, an eikonal approach to the shadowing corrections was proposed in the literature 
{Ll| . The starting point of these papers is the proof of the Glauber formula in QCD ||12| , 



which considers only the interaction of the fastest partons with the target. As in QCD we 
rather expect that all partons should interact with the target, a more general approach was 
proposed in ||. This approach was applied to the nuclear case in and to the nucleon 
case in |IO|jLl|| . Some of the main characteristics of this approach are its validity in a large 
kinematic region and that it provides a limit case of the GLR equation. In this paper the 
shadowing corrections will be estimated using the eikonal approach. In the next section we 
discuss this approach in more detail, pointing first our motivation. 

One of the most striking discoveries at HERA is the steep rise of the proton structure 
function F 2 (x, Q 2 ) with decreasing Bjorken x |Tj|. The behaviour of the structure function at 
small x is driven by the gluon through the process g — > qq. Therefore the gluon distribution 
is the observable that governs the physics of high energy processes in QCD. HERA shows 
that the deep inelastic structure function F%{x, Q 2 ) has a steep behaviour in the small- 
x region (10~ 2 > x > 10 -5 ), even for very small virtualities (Q 2 ~ lGeV 2 ). Indeed, 
considering F 2 oc x~ x at low x, the HERA data are consistent with a A that varies from 
0.15 at Q 2 = 0.85 GeV 2 to 0.4 at Q 2 = 20GeV 2 . This steep behaviour is well described 
in the framework of the DGLAP evolution equations with an appropriate choice of input 
distributions and the choice of the starting scale for the Q 2 evolution, by all groups doing 
global fits of the data [p}-|T6||. No other ingredient is needed to describe the experimental 
data. Using this result is it possible to conclude that there will be no significant shadowing 
contributions at HERA small- x regime? Here, we address this question. Clearly, the answer 
of this question is connected with the demonstration that the corrections due to shadowing 
corrections (SC) are negligible at least at HERA kinematic region. If it is not so the DGLAP 



approach is not better or worse than any other evolution mechanism developed to describe 
the experimental data. 

The SC for the F2(x,Q 2 ) proton structure function and xG(x,Q 2 ) gluon distribution 



were computed considering the eikonal approach in refs. ||10||ilf . In those works, the authors 
estimate the value of the SC which turn out to be essential in the gluon distribution but 
rather small in F 2 . Consequently, the F 2 data cannot determine if there is a new dynamical 
effect at HERA. Therefore, to obtain a more precise evidence of SC at HERA kinematic 
region, we must consider other observables directly dependent on the behaviour of the gluon 
distribution, and consequently, more sensitive to SC. Some of these observables are Fl, F 2 , 
J/ty production, diffractive leptoproduction of vector mesons and open charm production. 
In this paper we consider the observables Fl and F 2 C . The dominant charm production 
mechanism is the boson-gluon fusion. Therefore, the charm component F 2 of the structure 
function is a sensitive probe of the gluon distribution at small x. Furthermore, the violation 
of the Callan-Gross relation (Fl = F 2 — 2xF\ = 0) occurs when the quarks acquire transverse 
momenta from QCD radiation. Consequently, also the longitudinal structure function is a 
sensitive probe of the gluon distribution at small x. In this work we estimate the SC in Fl 
and F 2 in the approach proposed in [ID]. We calculate these observables using the Altarelli- 
Martinelli equation fl7| and the boson-gluon fusion cross section jL4J], respectively. The 
gluon distribution and the F 2 structure function necessary to these equations are calculated 
considering the eikonal approach. Actually, there are few data in these observables, but 
more measurements, with better statistics, will be available in the next years. Therefore it 
is very important the analysis of the behaviour of these observables considering SC. 

This paper is organized as follows. In Section II the eikonal approach and the shadowing 
corrections for the F 2 (x,Q 2 ) and xG(x,Q 2 ) are briefly presented. In Section III we obtain 
the formulae for the calculation of the Fl and F 2 and we apply our results to HERA data. 
Moreover, we compare our results with the predictions of DGLAP dynamics. Finally, in 
Section IV, we present our conclusions. 
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II. THE EIKONAL APPROACH IN PQCD 



The space-time picture of the Eikonal approach in the target rest frame can be viewed 
as the decay of the virtual gluon at high energy (small x) into a gluon-gluon pair long before 
the interaction with the target. The gg pair subsequently interacts with the target. In the 
small x region, where x « (R is the size of the target), the gg pair crosses the target 
with fixed transverse distance r t between the gluons. 

The cross section of the absorption of a gluon g* with virtuality Q 2 can be written as 



g.cf+nudeon^ 



Q2) = ( l dz \ d -^ \ ^W\Q 2 ,r t ,x,z)\ 2 a^ +nudeon (z,r?) , (1) 

JO J 7T J 7T 



where z is the fraction of energy carried by the gluon, 6 t is the impact parameter and \I/f 
is the wave function of the transverse polarized gluon in the virtual probe. Furthermore, 
a gg+nucieon^ z ^ r 2^ j g ^ e cross section of the interaction of the gg pair with the nucleon. In 
the leading log(l/x) approximation we can neglect the change of z during the interaction 
and describe the cross section a 99+nucieon^ z ^ r 2^ ag a f unc tj on Q f the variable x. 

Using the unitarity in the s-channel, the cross section (J aa+ nudeon ^ z ^ r ^ — a(x, rf) can be 
written in the form 

a(x, r 2 ) = 2 J d\ (1 - e -^(x,nA)^ ? ( 2 ) 

where Q is an arbitrary real function, which can be specified only in a more detailed theory 
or approach than the unitarity constraint. One of such specified model is the Eikonal model. 

The Eikonal model assumes that Q is small (Q << 1) and its bt dependence can be 
factorized as Q = QS(b t ), with the normalization J d 2 b t S(b t ) = 1. The factorization was 
proven for the DGLAP evolution equations and, therefore, all our further calculations will 
be valid for the DGLAP evolution equations in the region of small x or, in other words, in 
the double log approximation (DLA) of perturbative QCD (pQCD). The Eikonal approach 
is the assumption that Q = QS(b t ) in the whole kinematical region. 



In JT8| the authors demonstrate that is given by 
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{l = aff = ^^7r 2 r?xG(x, 7 ,) , (3) 

where xG(x, Q 2 ) is the gluon distribution. Therefore the behaviour of the cross section fl2|) in 
the small- a; region is determinated by the behaviour of the gluon distribution in this region. 
Considering the s-channel unitarity and the eikonal model, equation ([!]) can be written 

as 

(jg * + nuci e onu Q 2 ) = f 1 dz rfnf (Q*,r t ,x, z)\ 2 (1 - e^ ns W) , (4) 

JO J 71 J IT 

Using the relation a 9*+ nucleon (x^ Q*) = 47T ^ a xG(x, Q 2 ) and the expression of the wave \1> 9 * 
calculated in ||, the Glauber-Mueller formula for the gluon distribution is obtained as 

7T J Jx X' J-ky Tirf Jo 7T 

The use of the Gaussian parameterization for the nucleon profile function S(bt) = 

- h2 

simplifies the calculations. Consequently, doing the integral over b t , the master 
equation is obtained 

2R 2 /-i dx' r-gs d 2 r t 

where C is the Euler constant, E± is the exponential function, and the function Kc{x,r 2 ) = 
^rnr 2 xG(x, If equation (|6|) is expanded for small kq, the first term (Born term) will 
correspond to the usual DGLAP equation in the small x region, while the other terms will 
take into account the shadowing corrections. 

The master formula @ is correct in the double logarithmic approximation (DLA) ||10|| . 
As shown in |lfj the DLA does not work quite well in the accessible kinematic region 



xG(x,Q 2 ) = q- / ^ f^iC + Hnaix'y^ + E^aix'y,))} , (6) 
7r Jx r ii nrt 



(Q 2 > lGeV 2 and x > 10 -4 ). Consequently, a more realistic approach must be considered 
to calculate the gluon distribution xG. In (l(J the subtraction of the Born term of (||) and 
addition of the GRV parameterization was proposed. This procedure gives 



xG{x,Q 2 ) = xG master [eq. (6)] + xG (x, Q ) 

a s N c /-i f Q 2 dx' dQ' 2 x , g grv {x ^ g, 2) (7) 



7r Jx Jqi x' Q' 2 
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The above equation includes also xG GRV (x, Qq) as the initial condition for the gluon dis- 
tribution and gives xG GRV (x, Q 2 ) as the first term of the expansion with respect to Kq. 
Therefore, this equation is an attempt to include the full expression for the anomalous di- 
mension for the scattering off each nucleon, while the use of the DLA takes into account 
all SC. In || this procedure was applied to obtain the shadowing corrections to the nuclear 
gluon distribution and in |I(J to the nucleon gluon distribution. 

The equation (|) is not a non-linear equation type the GLR equation 0; it is the analogue 
of the Glauber formula, which provides the possibility to calculate the SC using the solution 
of the DGLAP evolution equation. In this paper, as shown in eq. ([?]), we use the GRV 
parameterization as a solution of the DGLAP evolution equation. It describes all available 



experimental data quite well |U|. It should be also stressed here that we disregard how 
much of SC has been taken into account in this parameterization in the form of the initials 
distributions. 

A similar approach can be used to obtain the SC for the deep inelastic structure function 
F 2 (x,Q 2 ). This was made in the references The main result of these references is 



that the structure function F 2 can be written, in the eikonal approach, as 

F 2 (x, Q 2 ) = ^ E 4 ^ 7? / *W ~ e^***)} , (8) 



where Vt qq = ^QS(b t ) (see eq. @). 

Following the same steps used in the case of the gluon distribution, the proton structure 
function can be written as 

OC>2 r l J2 

F 2 (x } Q 2 ) = E 4 ?° -^{C + ln(K q ( X y t )) + E 1 ( Kq (x,r 2 ))} , (9) 

u,d,s Q2 * 

where n q = |^g- Similarly as made in the case of the gluon distribution, to obtain a more 
realistic approach the Born term must be subtracted and the GRV parameterization must 
be summed. Therefore the proton structure function is given by 

F 2 (x,Q 2 ) = F 2 (x,Q 2 )[eq. (9)] - F 2 (x, Q 2 )[Born] + F 2 (x, Q 2 )[GRV] , (10) 
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where F 2 (x, Q 2 )[Born] is the first term in the expansion in n q of the equation (^|), and 

F 2 (x, Q 2 ) [GRV] = £ e\ [xq(x, Q 2 ) + xq(x, Q 2 )} + F 2 c (x, Q 2 ) (11) 

u,d,s 

is calculated using the GRV parameterization. 

In [[HJ the shadowing corrections for the F2(x,Q 2 ) proton structure function and 
xG(x, Q 2 ) gluon distribution were computed considering the eikonal approach, and the SC 
which turn out to be essential in the gluon distribution are rather small in F 2 . 

The result above is demonstrated in figure [l], where we present the ratios 

xG GM (x,Q*) 

1 xG GRV (x,Q 2 ) 1 ' 

and 

R2 ~ F GRV (x,Q*) ' (13) 

as a function of variable ln(-) for different values of Q 2 . In this case GM represents that 
the function (F 2 or xG) was obtained using the eikonal approach, and GRV represents that 
the function was obtained using the GRV parameterizations . We can see that in HERA 



kinematical region (3 < In- < 12) the behaviour of the gluon distribution is strongly 
modified by shadowing corrections. In |10| the authors demonstrated that the F 2 data can 



be described considering shadowing corrections, i. e. the F 2 data cannot determinate if there 
is a new dynamical effect at HERA. Consequently, observables directly dependent on the 
behaviour of the gluon distribution must be considered to discriminate the SC. From results 



obtained in [|1(| we conclude that the eikonal model provides a good description of the SC 
and can be taken as a correct approximation in the approach to the nucleon case for HERA 
kinematical region. In the next section we estimate the shadowing corrections to observables 
directly dependent on the behaviour of the gluon distribution. 

In this paper we estimate the shadowing corrections using an eikonal approach proposed 



in The eikonal approach gives sufficiently reliable results for the HERA kinematic 

region, however, it is not the most efficient way to calculate the shadowing corrections. In 
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T0|| , a generalized equation which takes into account the interaction of all partons in a 
parton cascade with the target was proposed. The main properties of generalized equation 
are: (i) the iterations of this equation coincide with the iteration of the Glauber-Mueller 
formula; (ii) its solution matches the solution of the DGLAP evolution equation in the 
double-logarithmic-approximation (DLA) limit of pQCD; (iii) it has the GLR equation as a 
limit, and (iv) contains the Glauber-Mueller formula. Therefore, the generalized equation is 
valid in a large kinematic region. In this paper we have used the Glauber-Mueller formula, 
since that in the HERA kinematic region the solutions of the generalized equation and 
of Glauber- Mueller formula are approximately identical [fTOfl. However, a more accurate 



approach, to a larger kinematic region than the HERA one, should consider as an input the 



solution of the generalized equation, in the future [19 



c 



III. F L AND F{ 

Our goal in this section is the determination of small- x shadowing corrections at Fl and 
F% using the approach proposed in jlTJ. These observables are directly dependent on the 
behaviour of the gluon distribution, and consequently, more sensitive to SC. 

Recently, the SC to the F2 slope was estimated |20]. This observable is also directly 



dependent on the behaviour of the gluon distribution ||21|| . It was shown that the contribution 
of SC is large in the F 2 slope and that the experimental data can be described incorporating 
the SC. This result is a strong motivation to estimate the SC at Fl and F 2 C . 

A. The longitudinal structure function 

The longitudinal structure function in deep inelastic scattering is one of the observables 
from which the gluon distribution can be unfolded. Its experimental determination is difficult 
since it usually requires cross sections measurements at different values of center of mass 
energy, implying a change of beam energies. A direct change of the beam energies at HERA 



has been widely discussed p2| . An alternative possibility is to apply the radiation of a hard 



photon by the incoming electron. Such hard radiation results into an effective reduction of 
the center of mass energy. Several studies on the use of such events to measure Fl have 



been carried out [2vJ. With these measurements, which in principle could be performed in 
the near future, it may be possible to explore the structure of Fl(x, Q 2 ) in the low x range. 
Current measurements of R(x,Q 2 ) = Fl/(F 2 — Fl) have been made by various fixed 



target lepton-hadron scattering experiments at higher x values ||24j| . At low x, the measure- 
ments of this observable are very scarce. Recently, the HI Collaboration has published the 
first Fl data at small x. To obtain the data, the structure function F 2 was parameterized 
taken only data for y < 0.35, where the contribution of Fl is small. The quantity 2/ = % 
describes in the rest frame of the proton, the energy transfer from the incoming to the out- 
going electron. This parameterization was evolved in Q 2 according to the DGLAP evolution 
equations. This provides predictions for the structure function F 2 in the high y region which 
allowed, by subtraction of the contribution of F 2 to the cross section, the determination of 
the longitudinal structure function. Therefore, in order to constrain the value of Fl in low x 
regime explored by HERA, the HI Collaboration extracted Fl assuming that the DGLAP 
evolution holds. Our point of view agrees with Thome's statement [|26||, pointing that the 



assumption of the correctness of the DGLAP evolution equations at small x implies that Fl 
is already isolated. Therefore, our comparisons with these data must be considered as an 
estimate of the SC compared with the results of the standard evolution equations. 

One can write the longitudinal structure function Fl in terms of the cross section for the 
absorption of longitudinally polarized photons as 

F L (x, Q 2 ) = * L (x, Q 2 ) - cr L (x, Q 2 ) (14) 

at small x. Longitudinal photons have zero helicity and can exist only virtually. In the 
Quark-Parton Model (QPM), helicity conservation of the electromagnetic vertex yields the 
Callan-Gross relation, Fl = 0, for scattering on quarks with spin 1/2. This does not hold 
when the quarks acquire transverse momenta from QCD radiation. Instead, QCD yields the 



Altarelli-Martinelli equation [17 
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F L (x, Q 2 ) = x 2 f %| F 2 (y, Q 2 ) + 4 £ e 2 (l - Q 2 )] , (15) 

An j x y 6 q y 

expliciting the dependence of Fl on the strong constant coupling and the gluon density. At 
small x the second term with the gluon distribution is the dominant one. Consequently, the 
expression (|T5| ) can be approximated reasonably by F L (x,Q 2 ) ~ 0.3 ^-xg(2.5x, Q 2 ) p7| . 
This equation demonstrates the close relation between the longitudinal structure function 
and the gluon distribution. Therefore, we expect that the longitudinal structure function 
to be sensitive to the shadowing corrections at HERA kinematic region. In this paper we 
calculate Fl using the Altarelli-Martinelli equation fll5]) . 

Considering the expressions (|7|), fllCf) and (|i~5|) we can estimate the SC for the F L structure 

function. In figure § we present the ratio 

_ F^(x,Q 2 ) 

~ F^ v (x,Q 2 ) ' (16) 
where F^ M (x, Q 2 ) indicates that the longitudinal structure function was obtained using 

the gluon distribution solution of the Glauber-Mueller formula, eq. ([5]), and F^ RV (x,Q 2 ) 

indicates that Fl was obtained using the GRV parameterization. We can see that the 

behaviour of Fl is strongly modified by shadowing corrections. The suppression increases 

with In(-) and is much bigger than for the Fi case. In the region of HERA data, 3 < 

ln( -) < 12, the shadowing corrections are not bigger than 44% (In(-) ~ 12). The SC are 

as big as 70% only at very small value of x (In(-) ~ 15), where we have no experimental 

data. 

In figure |3] we present the behaviour of the Fl structure function with (dashed curve) 
and without shadowing (solid curve) as a function of In(-) for different virtualities. We 



compare our results with the scarce HI data [25]. We can see that both curves agree with 



the data, however the experimental error are still very large. In figure |] we present our 
results for Fl and compare with all HI data. From the above results we can conclude 
that the eikonal model gives a good description of the longitudinal structure function and 
describes the experimental data. However, new data, with better statistics, could verify the 
shadowing corrections at HERA kinematic region. 
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Our conclusion of this section is that the longitudinal structure function Fl is a good 
observable to isolate the shadowing corrections at HERA. Consequently, we must stress the 
importance of measuring directly Fl at HERA as a probe of the dynamics of small x physics. 
We hope this result could also motivate the acquisition of new data in the next years. 



B. The charm component i 7 ^ of the structure function 

The problem of how to treat heavy quark contributions in the deep inelastic structure 



functions has been widely discussed, see for example p8| . It has been brought into focus 
recently by the very precise F%(x, Q 2 ) data from HERA [13]. Both the HI and ZEUS 



collaborations have measured the charm component F^ of the structure function at small x 
and have found it to be a large (approximately 25%) fraction of the total. This is in sharp 
contrast to what is found at large x, where typically F^jFi « 0(1O~ 2 ). 

For the treatment of the charm component of the structure function there are basically 
two different prescriptions for the charm production in the literature. The first one is 
advocated in [25| where the charm quark is treated as a heavy quark and its contribution is 



given by fixed-order perturbation theory. This involves the computation of the boson-gluon 
fusion process. In the other approach |30| the charm is treated similarly to a massless quark 
and its contribution is described by a parton density in a hadron. Here our goal is to obtain 
the shadowing corrections to the charm structure function in the eikonal approach, without 
considering in detail the question of how to treat heavy quark contributions ||28|| . However, 
some comments are important. We consider the charm production via boson-gluon fusion, 
where the charm is treated as a heavy quark and not a parton. This scheme is usually called 
Fixed Flavour Number Scheme (FFNS). In this scheme, by definition, only light partons 
(e.g. u, d, s and g) are included in the initial state for charm production: the number 
of parton flavours n/ is kept at a fixed value regardless the energy scales involved. The 
boson-gluon fusion gives the correct description of _F 2 C for Q 2 < Am 2 and should remain 
a reasonable approximation to F% for Q 2 > Am 2 . However, the boson-gluon fusion model 
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will inevitably break down at larger Q 2 values as the charm can no longer be treated as 
a non-partonic heavy object, and begins to evolve more like the lighter components of the 
quark sea. Therefore, our estimates in this scheme should be considered with caution in the 
region of large Q 2 . 

A cc pair can be created by boson-gluon fusion when the squared invariant mass of the 
hadronic final state W 2 > 4m 2 . Since W 2 = - + Mjy-, where is the nucleon mass, 
the charm production can occur well below the Q 2 threshold, Q 2 ~ Am 2 , at small x. The 
charm contribution to the proton structure function, in leading order (LO), is given by |TJ[ 

iF«z,QW) = 2e^> £^q 2 (^) 9 (y,^) , (17) 

where a = 1 + ^f- and the factorization scale \i! is assumed /i' 2 = Am 2 . G c g 2 is the coefficient 
function given by 

<W*, q! ) = 5 {(-- 2 + (i - -i 2 + »d " & > q# " * V 1 '" T=fi 

Am 2 

+ (3[-l + 8z(l-z)-z(l-z)-^]} , (18) 

where (3 = 1 — qt^^ is the velocity of one of the charm quarks in the boson-gluon center- 
of-mass frame. Therefore, in leading order, 0(a s ), _F 2 C is directly sensitive only to the gluon 
density via the well-known Bethe-Heitler process 7*g ~* c ~c- The dominant uncertainty in 
the QCD calculations arises from the uncertainty in the charm quark mass. In this paper 
we assume m c = 1.5 GeV. 

Considering the expressions (0) and ([[7]) we can estimate the shadowing corrections for 
the F% structure function. In figure |5] we present the ratio 



F^(x,Q 2 ) 
2 FZ GRV (x,Q 2 ^ ' 



We can see that the behaviour of the F£ is strongly modified by the SC. The suppression 
due to shadowing corrections increases with In(^) and is much bigger than for the F 2 case. 
In the region of HERA data, 3 < ln{~) < 12, the SC are not bigger than 40% {ln{\) « 12). 
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The SC are bigger than 62% only at very small value of x (In(-) ~ 15), where we have no 
experimental data. 

In figure |6| we present the behaviour of the _F 2 C structure function with (dashed curve) and 
without shadowing (solid curve) for different virtualities. We compare our results with the 
HI data [|3T]. Our predictions agree with the recent HI data. When compared with the GRV 
parameterization the shadowing corrections to _F 2 C present a suppression which increases with 
/n(^). Therefore, new data, with better statistics, could isolate the shadowing corrections 
at HERA kinematic region. 

As the shadowing corrections are distinct for the observables F2 and F% we can estimate 
the modifications in the charm contribution to the F 2 at small x by SC. In figure [7| we 
present our results for the ratio 

* = ffefi ■ (20) 

We can see that the behaviour of this ratio is strongly modified by the shadowing corrections. 
Therefore, we expect that the charm contribution do not increase at small values of x when 
compared with the predictions obtained using the usual parameterizations. Furthermore, 
we expect a strong modification in the charm production at small x. It is one of the main 
conclusions of this paper. This strong modification in the charm production shall occur in 
other associated observables, for instance in J/ty production. The shadowing corrections 



to the diffractive leptoproduction of vector meson J/ty was calculated in p2| . The authors 
obtained that the SC are large, but their results are affected by the uncertainty in the J/ty 
wavefunction. Consequently, although the cross section of diffractive leptoproduction of J/ty 
is proportional to the square of the gluon distribution, the accurate discrimination of SC in 
this process is still unlikely. 

Our conclusion of this section is that the charm component F£ is also a good observable 
to isolate the shadowing corrections. Therefore, new data, with better statistics will be very 
important in the determination of dynamics at HERA, and further applications. 
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IV. CONCLUSIONS 



The pQCD has furnished a remarkably successful framework to interpret a wide range 
of high energy lepton-lepton, lepton-hadron and hadron-hadron processes. Through global 
analysis of these processes, detailed information on the parton structure of hadrons, es- 
pecially the nucleon, has been obtained. The existing global analysis have been performed 
using the standard DGLAP evolution equations. However, in the small x region the DGLAP 
evolution equations are expected to breakdown, since new effects must occur in this kine- 
matical region. One of these effects is the shadowing phenomenon. 

The recent F 2 data are well described in the framework of the DGLAP evolution equa- 
tions with an appropriate choice of input distributions and the choice of the starting scale 
for the Q 2 evolution. Although no other ingredient has been needed to describe the data, 
we cannot conclude that the shadowing corrections are negligible at HERA kinematic re- 
gion only from the analysis of F 2 data. The proton structure function is inclusive on the 
behaviour of the gluon distribution. In this paper we estimate the shadowing corrections to 
Fl and F 2 °- The behaviour of these observables is directly dependent on the behaviour of 
the gluon distribution and, therefore, strongly sensitive to the shadowing corrections. We 
shown that the SC to Fl and F 2 are much bigger than for the F 2 case. Our results are in 
accord with the recent few data. New data with better statistics are strongly welcome. 

One of the main conclusions of this paper is the strong modification of charm component 
of the structure function. As the SC are different in the observables F 2 and F 2 the charm 
contribution do not increase at small x as predicted by the usual parameterizations. 

In this paper we estimate the shadowing corrections using an eikonal approach proposed 



in [ 10[ . The eikonal approach gives a sufficiently reliable result for the HERA kinematic 
region, however, a more accurate approach, to a larger kinematic region than the HERA 
kinematic region, should consider as an input the solution of the generalized equation pro- 
posed in |TIJ. In a future publication we will estimate the behaviour of some observables 
using the solution of the generalized equation. 
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The determination of the dynamics at small x is fundamental to estimate the cross 
sections of the processes which will be studied in the future accelerators. We expect that 
this paper could motivate a more accurate determination of Fl and F 2 C in the next years, 
since the behaviour of these observables may explicitate the dynamics at small x regime. 
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FIG. 1. (a) The ratio R x ■ 
ln{^-) for different virtualities. 



xG aM 
xG^V 



and (b) the ratio i?2 = F anv as a function of the variable 
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FIG. 2. The ratio Rl = f orv as a function of the variable In(^) for different virtualities. 
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FIG. 3. The longitudinal structure function Fl for Q 2 = 8.5, 12, 15GeV 2 . The solid curve 
represents that the Fl was calculated using the GRV parameterization and the dashed curve that 
its was calculated using xG and F<i obtained in the eikonal model. Data from HI. 
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FIG. 5. The ratio R% = 2 grv as a function of the variable In(^) for different virtualities. 
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FIG. 6. The charm component of the structure function for different values of Q 2 . The solid 
curve represents F| calculated using the GRV parameterization and the dashed curve as calculated 
using xG obtained from the eikonal model. Data are from HI. 
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FIG. 7. The ratio Rp = p 2 - as a function of the variable In(-) for different virtualities. The 
solid curve represents the predictions obtained using the GRV parameterization and the dashed 
curve the predictions obtained in the eikonal model. 
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